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^ . We study numerically magnetic Gm{p) and electric Ge{p) gluon propagators 

^ ! and their dependence on the choice of the lattice gauge field Ax/j. in SU (2) gauge 

O I theory, especially, in the low-momentum limit. We find that two different Ax^ 

• definitions are equivalent up to a trivial renormalization of the propagator, at least, 

■ in the main approximation. 

'f'- '. 1 Introduction 

> ■ 

•i-H . 

^ ' Gauge variant Green functions are supposed to contain an important information 

■ about the large distance physics and mechanism(s) of confinement. For example, at 
zero temperature the expected infrared suppression of the gluon propagator in the 
Lorentz (or Landau) gauge is believed to be connected to the Gribov's confinement 
scenario 0, 0] (see also 0). In the high temperature (chromoplasma) phase the 
large distance and low-momentum dependence of the gluon propagator may provide 
valuable information about the electric and magnetic screening mechanisms |^. 

On the lattice the definition of gauge variant Green functions, e.g. gluon prop- 
agators, is not free from the ambiguity. Indeed, the lattice gluon propagator G/^n, is 
given by 



G^^ix -y)^ (Tr [Ax^^Ay,] ) , (1.1) 

where Axfi are gauge fields on the lattice and (...) means the statistical average 
with some gauge fixing condition. There is no unique definition of Ax^ on the 
lattice in terms of the link variables Ux^- This problem has been addressed recently 
in 1^, ^, ^ and some possible choices of Axfi have been discussed. 

* Work supported by the grant INTAS-00-00111 and RFBR grant 02-02-17308. 



1 



A similar ambiguity exists in the choice of the gauge which is the Lorentz gauge 
in our case. On the lattice there are different ways to choose the Lorentz gauge 
fixing condition, all of them being the same in the continuum limit. 

In this paper we study the space-space (magnetic) Guip) and time-time 
(electric) Ge{p) correlators in the high temperature plasma phase with a special 
attention paid to their low-momentum behaviour. The main point in our study is 
the dependence of the correlators on the choice (definition) of lattice gauge field Ax^i 
as well as on the specific lattice implementation of Lorentz gauge fixing condition. 

In what follows the periodic boundary conditions are used, the lattice size is 
V4 = A^4 X and (9^ is a backward lattice derivative. We use also IC^{p) = 

1 sin ^ and /C^(p) = J^^i^^fiip) , ci being a lattice spacing. 

2 Gauge fields and Lorentz gauge conditions 
2.1 Definitions of the gauge fields 



The standard Wilson action [jTO| with SU (2) gauge group is 



/5 = 4> (2-1) 



where g is a bare coupling constant and U^^ G SU{2) are link variables. Under 
gauge transformations VL field variables Ux^i transform as follows 

Ux^. ^ U% = QxUx^nU^ ; Qx G SU{2) . (2.2) 

A standard definition [|ri|] of the lattice gauge field Axfj. in terms of link variables 
Ux^, is 

■^- = 2^K-^-)' (2.3) 

However, this definition is not unique. For example, instead of field A^^ one can 
define gauge field A' 



^^•. = i^K-(f^i.)')' (2-4) 
or another gauge field v4^^, 

etc. On the other hand, the proof of the (naive) continuum limit existence of the 
lattice action S{U) as well as the construction of the lattice perturbation theory is 
based on the representation of the link variable Ux^ in the form 
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= exp^^iagA^^^ , \Axf,\ < — , (2.6) 

cig 

and the expansion in series in powers of coupling constant g. Therefore, from 
the point of view of the analytical (e.g., perturbative) studies, there is a natural 
definition of the lattice gauge field, Ax^i 

45 = Ax.iU.,) , (2.7) 

Ax^ being defined in eq. (|2.6|) . Given any group element U = Col + iac = exp{^9a} , 
one can easily find that 

- 2c , , 

9 — -j-ij- arccos Co , (2.8) 

unless Co 7^ —1 (the case Co = — 1 has measure zero and can be discarded). A 
formal expansion in powers of spacing a gives 

<^ = 4? + 0(a2), (2.9) 

so that both definitions are equivalent in the naive continuum limit. 
In the rest of this paper the lattice spacing a is chosen to be unity. 

2.2 Choice of the Lorentz gauge fixing condition 

In lattice calculations the usual choice of the Lorentz gauge condition is [ITT 



E MS = 0' (2-10) 



which is equivalent to finding an extremum of the functional Ffj (fi). 



4'^(^) = 77rE^Tr[/f;,, (2.11) 

with respect to gauge transformatins f^a, . In what follows this gauge is referred to 
as LGq. Equally, one can choose another form of the Lorentz gauge fixing condition 
(see also 0) which will be referred to as LG\ : 

EMS = 0. (2.12) 

Evidently, both LGq and LG\ are the same in the continuum limit. The LG\ gauge 
condition is equivalent to finding an extremum of the functional fIj\Q,) 



XjJ, 
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The field A^fj, = Axfj.{Uxfj.) shows the nonanalytic dependence on the hnk variable 

(2) 

Ux^- However, the functional Fj) (Q) still remains a continuous function of Q. 
Under the infinitesimal gauge transformation SQ^ = exp{iSuJx} the variation of 
the functional fIj\^1) is 

5Ff (5n) = -^^^5^," + ... , (2.14) 

a 

where 

V.S = EMSm(^-m) ■ (2-15) 



Eqs. ( p.l4|) , (p.l5|) define the numerical algorithm of gauge fixing. To maximize 

fIj\^1) one can choose the gauge transformation matrix fi^, = exp{iuJx}, 

ul = -h-^l- b>0, (2.16) 

successfully at all lattice sites x. After a certain number of gauge fixing sweeps a 
local maximum Fj^lx{U) of the functional fIj\Q) is reached with given accu- 
racy. The value of the parameter b in Eg. ( p. 161) should be tuned to optimize the 
convergence. 

2.3 Lorentz gauge and Gribov copies 

It is well-known that gauge fixing on the lattice as well as in the continuum is 
affected by the existence of Gribov copies [Q. If one repeatedly subjects a configu- 
ration {Uxfi} to a random gauge transformation as in Eq. (|2.2|) and subsequently 
applies to it the Lorentz gauge fixing procedure, one usually obtains Gribov (or 
gauge) copies with different values of FmaxiU). It is a long-standing believe 
that the "true" gauge copy corresponds to the absolute maximum of Fu{VL). 

In our case the question of interest is the dependence of lattice quantities in the 
Lorentz gauge (s) on the choice of the gauge copy. This dependence has been well 



studied in a number of papers (see, e.g. |T^, |T4|, [T5|, |1^, [T^l)- It has been found that 



for different choices of gauge copies variation of the gluon propagator is small and 
comparable with statistical error. 

Our study confirms this statement. For every equilibrium configuration {U^^} 
we performed 10 random gauge transformations and then applied a gauge fixing 
procedure. For the calculation of the gluon propagator we have used either the 
first gauge copy of the equilibrium configuration {U^^i} or the "best" one, i.e. the 
copy with the maximal value of FmaxiU). The difference between the (averaged) 
propagators appeared to be very small and negligible in comparison with statistical 

^It is worth noting that in some of the above mentioned papers many more gauge copies have 
been used. 
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errors. Therefore, we conclude that there is no in fact Gribov copy problem, at 
least, for the SU{2) lattice gluon propagator in the Lorentz gauge. 

Note that this situation is quite different from what one observes in the case of 
[/(I) lattice gauge theory. Indeed, the id photon f/(l) propagator in the Lorentz 
gauge exhibits an extremely strong dependence on gauge copies JTBI, IT^ . 



3 Numerical results 

We calculated magnetic Gm{p) and electric Ge{p) correlators, 

Gm{p) = l{Gn{p) + GM) ; Ge{p) = GM , (3.1) 
where G^y{p) is given by 

GAp) = ^ (Tr (A(p)A(-p))) . (3.2) 
Field A^{p) is a Fourier transform 

^A.(p)=Ee-^''"^"^-A.;., (3.3) 



where Ax^j. are defined in Eqs. ( |2.3| ) and ( p. 61) with unit spacing. Momenta p 
have been chosen to be directed along the third axis, i.e. p = (0,0,^3,0) with 
P3 = and ns = 0, 1, .... 

Most of our calculations have been performed on 4 x lattices with Ng = 
16; 24; 32; 40 and on the x 24^ lattices with = Q; 8 . The following 
values of (3 have been chosen: (3 = 2.512; 2.645 and 2.74. The choice of (3 = 2.512 
corresponds to the temperature T = 2Tc on the lattice with A^4 = 4 as well as the 
choice of (3 = 2.645 on the lattice with N4 = 6 and /3 = 2.74 on the lattice with 

= 8 i3 1 



As it has been already mentioned, at zero temperature one expects the infrared 
suppression of the gluon correlator, i.e. the nonperturbatively calculated gluon cor- 
relator is expected to be less singular than the perturbative one. It has been shown 
that Lorentz (and Coulomb) gauge fixed gluon correlators on the infinite lattice 
vanish in the zero momentum limit due to the proximity of the Gribov horizon 
(This has been confirmed recently in the framework of Dyson-Schwinger equation 
formalism 0). Recently this study has been extended to the finite temperature 
case |]21|. It has been shown that in the infinite volume limit magnetic correlation 



functions vanish in the infrared unlike the electric correlation functions . 

In Figure |l| one can see the momentum dependence of the magnetic correlator 
G^'*(p) for different lattice sizes. At larger values of the momentum (/C^^O.2) 
the finite volume dependence is practically absent. On the other hand, at values 



■'■We are grateful to F. Karsch and J. Engels for providing us the /3-value for N4 = 6 
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/C^ ~ this dependence is rather strong. For our largest lattice with A^^ = 40 
the momentum dependence acquires a local maximum at some nonzero value of p. 
A similar structure in the high temperature phase has been observed recently in 
||2^ (compare with the low temperature case 0, In Figure |^ we show 

the dependence of the zero momentum correlator G^^j{0) on A^^ at A^4 = 4 and 
f] = 2.512. One can see a monotonic descrease of Gm{0) with increase of the lattice 
size Ng. However, one needs much larger lattices to judge on the infinite volume 
limit. 

The momentum dependence of the electric correlator G^^^ (p) shown in Figure ^ 
differs strongly from that of the magnetic correlator, in agreement with [^. All 



curves are monotonic and there is no any local maximum at nonzero momentum 
p . Note also that at small momenta (/C^ ~ 0) the finite size dependence of G^^\p) 
is rather weak, in contrast to G^lj{p) case. 

To see the dependence of the correlators on the choice of the lattice gauge fields 
we defined ratios Rm and Re '■ 

These ratios are shown in Figures ^,b for different lattices at /? = 2.512 , LGq and 
LGi. Up to "second order" corrections (not visible in these Figures) both ratios 
are equal and momentum independent: 

Rm{p; g'; T) ~ Re{p-, g"-, T) :^ C{g^; T) . (3.5) 

A more detailed structure is shown in Figure § for the ratio of magnetic correlators 
under LGq and LGi at fixed temperature T = 2Tc and different /?'s. One can see 
rather weak (~ 1%) dependence on j3 and even weaker dependence on the choice of 
the gauge condition (i.e. LGi vs. LGq). Another observation which can be made in 
Figure ^ is a - very weak - momentum dependence of the ratios. At the moment, 
we cannot exclude that this dependence is a numerical artifact. This point needs 
more detailed study. 



4 Conclusions 

We studied numerically the transverse magnetic Gm{p) and electric Ge{p) gluon 
propagators in the high temperature plasma phase (T > T^) in the SU{2) lattice 
gauge theory. 

The transverse magnetic propagator Gm{p) exhibits very strong volume de- 
pendence in the infrared region. The value Gm{0) tends to decrease with volume 
increase in agreement with the prediction by Zahed and Zwanziger It is worth- 



while to note that one needs rather large lattices to study this effect in details. 
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The momentum dependence of the electric correlator is very much different 
from that for the magnetic correlator. There is no any local maximum at nonzero 
momentum p and in the infrared region the finite size dependence of Ge{p) is 
rather weak, in contrast to Gm{p) case. 

We studied the dependence of the propagators on the choice of the lattice gauge 
field Axfj., one of them being standard and another one being "natural" from the 
viewpoint of the perturbation theory, as well as their dependence on the choice of the 
lattice Lorentz gauge. We found that two different Ax^ definitions are equivalent up 
to a trivial renormalization of the propagator, at least, in the main approximation. 
Most probably, this proportionality factor can be explained (at least, partially) by 
the tadpole renormalization (see also the discussion in p5|). 
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Figure 1: (p) at (3 — 2.512 for LGq. The line is to guide the eye. 
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Figure 2: The dependence of G^''(0) on Ng aX (5 = 2.512 for LG^. The hne is to 
guide the eye. 
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Figure 3: G%\p) at /3 = 2.512 for LGq. The line is to guide the eye. 
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Figure 4: Ratios of magnetic (a) and electric (b) correlators at /5 = 2.512 for LGq. 
Lines are to guide the eye. 
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Figure 5: Ratios of magnetic correlators at /3 = 2.512 
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